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ABSTRACT 

I consider the possibility that Ultra High Energy Cosmic Rays are accelerated in 
Gamma Ray Bursts located in the Galactic corona, thus circumventing the problem 
raised by Greisen-Zatsepin-Kuz'min cutoff. The acceleration of UHECRs could occur 
in the pulsars which, in the coronal GRB model, produce them: the same parameters 
that permit fitting GRBs' observations in the model of Podsiadlowski, Rees and Rud- 
erman (1995) lead to an estimate of the highest achievable energies corresponding to 
that of the Bird et al. (1994) event, and to very low luminosities in cosmic rays. I show 
that, if the observations of Milgrom and Usov (1995a) are confirmed, the extragalactic 
GRBs' model for the acceleration of UHECRs is untenable, but the same constraint 
does not apply to the coronal model. Also, I show that the efficiency of particle accel- 
eration needs be much smaller (and less demanding) than in cosmological models of 
GRBs. Uncertainties remain about the ensuing cosmic ray spectral distribution. I also 
briefly discuss observational strategies to distinguish between the two possibilities. 
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1 INTRODUCTION 



The recent discovery of the two highest energy CRs ever 
(Bird et al, 1994, Yoshida et al., 1995) has produced a 
renewal of interest in the acceleration of Ultra High En- 
ergy Cosmic Rays (UHECRs). This is due to the well- 
known Greisen-Zatsepin-Kuz'min effect: because of photo- 
pion losses, a proton with E = 10^^ eV at the source needs 
only cross a distance ~ 30 Mpc (Protheroe and Johnson 
1995) to be slowed down to S = 3 x 10^° eV. Thus it is 
necessary to identify acceleration sites which are reasonably 
close to the Earth. 

It has been recently pointed out (Milgrom and Usov 
1995a, Waxman 1995a, Vietri 1995) that UHECRs may be 
accelerated in cosmological Gamma Ray Bursts, GRBs. It 
was shown (Vietri 1995) that CRs'energies as high as the 
largest ever observed {E « 3x 10^° eV, Bird et al, 1994) can 
be achieved in the very short durations (« 1 sec) of GRBs 
by means of very efficient first-order Fermi-Bell accelera- 
tion in just two cycles. However, as it is well-known (Lamb 
1995, Paczyriski 1995), the nature and location of GRBs is 
in dispute, the most likely alternative to the extragalactic 
model being a distribution of neutron stars in an extended 
Galactic corona, at Galactocentric distances « 100 kpc, and 
it may be interesting to consider whether UHECRs can be 
accelerated in coronal GRBs. 



This idea seems to subvert the traditional view (Cocconi 
1956) that UHECRs are of extragalactic origin, the main ev- 
idence for this lying in the change of chemical composition 
and spectral slope at the 'ankle' (_E = 3 x 10^* eV). However, 
the objects postulated to give rise to coronal GRBs form a 
unique population, not just because their flux and angular 
distributions are sharply at odds with those of all known 
Galactic populations and identical to those of all known ex- 
tragalactic sources, but also because their Galactocentric 
distances exceed by an order of magnitude those of known 
Galactic objects. Thus, the hiatus that separates Galactic 
and extragalactic UHECRs' sources exists in this model as 
well. The classical way to try to establish the extragalac- 
tic nature of the sources giving rise to UHECRs has been 
to seek the GZK-cutoff in the CRs' flux observed at the 
Earth. At present, however, the evidence is lacking statisti- 
cal significance (Sigl et al., 1995), and does not rule out the 
coronal model. 

It seems thus worthwhile to consider whether UHECRs 
can be originated in these nearer models for GRBs. I shall 
point out below three reasons why this seems attractive. A 
discussion will follow in Section 5. 
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2 THE ACCELERATION OF UHECRS IN 
CORONAL GRBS 

The essential feature of the extragalactic GRB model that 
is preserved by the Coronal model is that GRBs are lo- 
cated inside the Greisen-Zatsepin-Kuz'min (GZK from now 
on) sphere. The traditional acceleration sites for UHECRs, 
Fanaroff-Riley Class II radio galax;ies (Rachen and Bier- 
mann 1993) certainly provide an energetically attractive 
source of UHECRs, despite the large uncertainties in the es- 
timates of the proton flux at the source. However, the nearest 
such galaxy, Pictoris A, lies at dpA ~ 100 Mpc. The GZK- 
radius ioi E ^ 2 x 10^" eV is Rgzk = 20 Mpc (Protheroe 
and Johnson 1995), so that the total flux emitted by Pictoris 
A is damped by the factor exp{dpA/RGZK) « 100. It will 
be shown later that 1 extragalactic GRB is expected within 
a GZK-sphcrc every 10'^ years. If this emits equal amounts 
of energy in 7-ray photons and CRs (not just UHECRs), 
then the expected ratio of fluxes at the Earth is 

« 0.1 expidp a/ Rgzk) « 10 . (1) 

JPicA 

The super-GZK flux by FRII galaxies is oven more strongly 
dominated by coronal GRBs: all coronal GRBs are located 
inside Rgzk, as opposed to just a fraction {RgzkHq/c)^ « 
10~^, giving a flux of super-GZK CRs higher by a factor 
c/HqRgzk ~ 10^. Also, coronal GRBs, exactly like their ex- 
tragalactic counterparts, are largely super-Eddington, and 
hyperrelativistic phenomena like beaming naturally arise 
around them, thus mirroring the discussion in Vietri (1995) 
that made GRBs attractive as potential sources of UHECRs. 

The mechanism for the acceleration of UHECRs in cos- 
mological GRBs (Vietri 1995) surely docs not work for the 
coronal model, because the model of an extragalactic GRB 
as due to the prompt release of a shock's whole energy im- 
mediately after formation, by self-synchro-Compton of rel- 
ativistic electrons (Meszaros and Rees 1994) requires high 
ISM densities to achieve sufficient ofBcicncics, n ~ 1 cm^' 
(see also Begelman, Meszaros, Rees 1993). Such high bary- 
onic densities are inconceivable at the distances (« 100 kpc, 
Podsiadlowski, Rees and Rudorman 1995, PRR from now 
on) currently postulated for the coronal scenario, where the 
dark matter density is pdm ~ 3 x 10~^mH cm~'^, and thus 
most likely n ^ lO"** cm~^. Models for coronal GRBs can- 
not simply be those concocted when GRBs were thought to 
lie at « 1 kpc from us, because the greater distance scale im- 
plies release of ~ 10* times more energy. Recently, however, 
PRR have shown that a reasonable GRB-generation mech- 
anism can be identified in the stress-release episodes of the 
crustal magnetic field, provided B « 10^® G. This suggests 
that cosmic rays may be accelerated in pulsars' magneto- 
spheres. Sigl, Schramm and Bhattacharjee (1994) give as 
the highest cosmic ray energy from a pulsar 

S^a. = 2 X 10^° eV (yoJ^) • (2) 

The above formula was discussed in this context also by Mil- 
grom and Usov (1995a); the novel point is the 'coincidence' 
that the same magnetic field is necessary to explain both the 
GRB and the UHECR phenomena. While no more specific 
predictions can be made because of the lack of a detailed 
coronal GRB model, I find this coincidence encouraging. 



3 COINCIDENCES BETWEEN UHECRS AND 
GRBS 

Another reason why the coronal gamma-ray-burst hypoth- 
esis is attractive comes about when we consider the impli- 
cations of the work by Milgrom and Usov (1995a). They 
discovered that the two super-GZK events (Bird et a/. 1994, 
Yoshida et al, 1995) were positionally coincident, within 
their largish error boxes, with two strong GRBs which 
preceded them hy ^ 1 yr; on the basis on this associa- 
tion they proposed that UHECRs are generated in GRBs. 
If one believes in this association, the following argument 
shows that GRBs are unlikely to be extragalactic. Let us 
compute the expected rate of GRBs resulting in super- 
GZK events at the Earth. Given that the rate of GRBs is 
(Paczyriski, 1993) h — 30 Gpc^'^ yi~^^ , and that a Greisen- 
Zatscpin Kuz'min sphere with radius Rgzk = 20 Mpc, 
for £; = 2 - 3 X 10^° eV, (Protheroe and Johnson 1995) 
has volume Vgzk = 3 x lO'* Mpc^, I find that a GZK- 
sphcrc has a rate of 1 GRB every lO"* yr. However, in the 
« 10 yr of combined observational time of the Fly's Eye and 
AGASA experiments, two such events have been already ob- 
served. This occurs, in the above model, with probability 
P2 = (lOg/10^)^ « 10"^ where q is the averaged fraction 
of sky coverage, which has been taken as g ^ 0.3. Stated 
another way, for these low values of Rgzk and of the time- 
delay, most of the time (always but just once in 10^ yr) we 
should not be able to observe any super-GZK CR. It follows 
that, if the proposed observational connection and time de- 
lay (^ 1 yr) were to bo confirmed, the sources that produce 
UHECRs could not be extragalactic GRBs. 

The above observation is of course unlikely even when 
viewed from the other angle, that of GRBs. Milgrom and 
Usov (1995a) had at their disposal ~ 3 j/r of the BATSE cat- 
alog. This implies then that the probability of seeing 2 GRBs 
located within Rgzk is also given by (3/10^)^ = 10~^, for- 
tuitously equal to P2 above. Another way to look at the 
same problem is to compute the total energy released by 
the two GRBs that Milgrom and Usov (1995a) associated 
with the two highest energy CRs. They have fluences of 

4 X 10~^ erg cm~^ and 3 x 10~* erg cm~'^, which assuming 
a distance = Rgzk = 20 Mpc, correspond to a total energy 
release of lO*** — 10*'' erg. This is low when compared with 
the average energy released by GRBs, 4 x lO^'^ erg (Piran 
1992). While it is certainly possible that the GRBs' luminos- 
ity function is broad, still the smallncss of the total energy 
released computed thusly is in keeping with the argument 
developed above. 

It should be noticed that the previous argument is based 
upon the smallness of Rgzk = 20 Mpc. Protheroe and 
Johnson (1995) compare their results with several previous 
computations, and it is apparent from their Fig. 4 that theirs 
is the largest value in the energy range of interest here, some 
authors having obtained values as low as Rgzk ~ 9 Mpc. 

Independently of the actual time-delay, it seems very 
unlikely to me that any GRB in the BATSE catalog can 
come from within the GZK-sphere. One way out of this 
predicament is, of course, if the time-delay is actually much 
longer than proposed by Milgrom and Usov (1995a); in this 
case there need be no GRB from inside the GZK sphere in 
the BATSE catalog, and, if it is assumed that CRs from the 
same source are spread out over a time comparable to the 
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time-delay, then perhaps a hundred distinct extragalactic 
GRBs as sources of UHECRs should be seen at all times (see 
later on). More data should settle this issue. In the following 
I shall Eissume the connection and the time— delay ~ 1 yr to 
be correct. In this case the above predicament would be 
relieved if coronal GRBs were to produce super-GZK CRs, 
first because photopion and photoelectron losses within the 
Galaxy are entirely negligible (and thus several GRBs from 
within the UHECR error box can be the putative fathers of 
the super-GZK events), and second because there is no need 
that the objects giving rise to GRBs produce UHECRs only 
during fireballs: they could instead produce a steady fiux 
of CRs. Under the coronal hypothesis, the time delays of 
« 1 j/r can easily be accommodated (Milgrom and Usov 
1995a). 



4 THE EFFICIENCY 

Another reason why a Coronal GRB origin for UHECRs 

is attractive is due to the fact that it was quickly realized 
(Vietri 1995, Waxman 1995a, Milgrom and Usov 1995b) 
that accounting for all UHECRs observed at Earth re- 
quires that each extragalctic GRB releases approximately 
equal amounts of energy in 7 ray photons and in UHE- 
CRs, i.e., very high efficiency. In order to see how seri- 
ous this efficiency problem is, it is convenient to compare 
the expected production by GRBs of UHECRs in the range 
10" eV < E < 10^° eV with that deduced from observa- 
tions (Waxman 1995b), eohs = 5 x 10~^^ erg cm~^. The 
observed rate of GRBs is h = 30 Gpc^''' yr~^ (Paczynski, 
1993). The average GRB total energy output in photons E-^ 
is 4 X 10®^ erg (Piran 1992). This yields the total energy 
release rate in 7-ray photons. To estimate the release rate 
of energy in UHECRs I proceed as follows. The fraction q 
of the directed kinetic energy that shocks convert into CRs, 
compared to that converted into thermal motions (which in 
GRBs is promptly dissipated into 7-ray photons) is very 
uncertain, and is variously estimated (Draine and McKee 
1993) between 0.03 and 0.5. I shall take the larger value. The 
spectrum emitted by a cosmological source is derived from 
observations (Waxman 1995b) as oc E~^"^ . I shall consider 
a harder spectrum, oc E~'^ , so as to continue to overesti- 
mate the production of UHECRs by GRBs. Then the total 
fraction of energy released in cosmic rays, and channeled in 
the range discussed by Waxman, 10^® eV < E < 10^" eV , 
is p = hi{Ejnaxl Emin)- I shall conservatively take 

Err,ax = 10^° eV , and S^i™ = 10^^ eV . This last value 
comes from this argument. Vietri (1995) took as would 
be CRs the extreme Boltzmann tail of the just-shocked 
ISM protons, for which 7 w 10^ — 10* in the shell frame. 
In the lab frame these protons would appear as CRs with 
E ~ lO'^* — lO'^^ eV , corresponding to the range in 7. To- 
gether, these two factors imply a relative efficiency of UHE- 
CRs to 7-ray photons of ?7 = Euhecr/Ej = qp k, 1/10. 
The comparison with observations yields 



^obs ^obs 

despite my attempts at mgiximizing the contribution of 
GRBs. Thus it seems likely that, in order to reproduce ob- 
servations, GRBs must overproduce UHECRs with respect 



to conventional models, by a factor 10 — 100, i. e., ~ 1 — 10. 
It should be emphasized that large efficiencies are not impos- 
sible: the above estimate is very uncertain, and it refers to 
steady-state, newtonian shocks because no analogous com- 
putations are known to me in the time-dependent, relativis- 
tic regime. However, in the face of the daunting task of rais- 
ing the efficiency by about two orders of magnitude, it seems 
worthwhile to consider the alternative hypothesis that GRBs 
originate in the Galactic corona. 

The efficiency requirement becomes immediately less 
stringent in the coronal model. In fact, if UHECRs are gen- 
erated by cosmological GRBs, only those originating within 
RgzK; the Grcizen Zatsepin Kuz'min radius, of the Milky 
Way can reach us. This is because UHECRs with E ^ & few 
times 10^'' eV loose energy by photopion and photoelectron 
production off CMBR photons. These nearby GRBs account 
for a fraction / « Rgzk/{c/Ho) of the whole fiux at Earth. 
On the other hand, if GRBs are located in the Galactic 
Corona, all GRBs generate UHECRs that reach the Earth. 
The coronal GRBs' energy release (Egrb ~ 10*^ erg) is de- 
termined so that the flux of photons at Earth is equal to that 
in the cosmological model; thus the coronal model with the 
same efficiency factor rj as the extragalactic model produces 
an UHECRs' flux at Earth higher by the factor 1/f with 
respect to the extragalactic model. With Rqzk = 20 Mpc 
(Protheroe and Johnson 1995), for S w 3 x 10^° eV , the 
overproduction is 1// « 10^. This implies that, in the coro- 
nal model, the fitting of the observed UHECRs' flux can 
be achieved with an efficiency r; reduced by the factor / to 
7? « 0.03. 



5 DISCUSSION 

The low efficiency rj « 0.03 discussed above is rather re- 
warding in the case in which UHECRs are generated dur- 
ing GRBs. However, another, more effective way for the low 
value of / to ease our luminosity quandaries occurs if the 
UHECR-production is steady, within the PRR model. In 
fact, PRR postulate that every pulsar remains active for up 
to lO^'' yr, producing ~ lO" stress release episodes within 
its lifetime. This means a GRB every lO'* yr, which equals a 
time-averaged GRB-luminosity of w 3 x 10^® erg . With 
rj « 0.03, this leads to a continuous UHECR-luminosity of 
LuHECR ~ 10^* erg s~^. 

However, let me remark that, if the release of the UHE- 
CRs were coincident with the event leading to the GRB, the 
acceleration of UHECRs could use as an energy source the 
very same one of the GRB, and it would be energetically in- 
significant since rj ~ 0.03. This would also leave room for the 
acceleration of several other decades of CRs' energy beyond 
the one considered, 10^^ eV < E < 10^° eV, even assummg 
a softer spectrum. It is more difficult to identify an energy 
source in the case of continuous acceleration of CRs: with a 
field B ~ 10^^ G, the rotational kinetic energy of the pulsar 
must have been exhausted very early indeed. 

Another observational feature, the dominance of pro- 
tons with respect to heavier nuclei beyond the ankle ( 
£ = 3 X 10^* eV, Bird et at, 1994, 1995), which is often 
cited as evidence for the extragalactic origin of UHECRs, 
can be accommodated easily within this model. Iron nuclei 
with E fs 10^° eV have gyroradii rz, « 5 kpc in the Galactic 
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magnetic field, so that they are essentiaUy confined around 
the pulsars producing them. All diffusion processes tend to 
push them away from us, and into the IGM, thus making 
them unobservable to us. Protons of comparable energy have 
much larger gyroradii (« 100 kpc), which allow them to pen- 
etrate into the inner Galaxy. 

The most significant weakness of the coronal model is 
that the expected spectrum is not necessarily close to E~'^, 
as it automatically is every time OR acceleration at shocks is 
invoked. This point is clearly in need of further investigation. 

Neglecting the thornier question of establishing a con- 
nection between UHECRs and GRBs, I discuss now a com- 
parison between the cxtragalactic and coronal models for 
GRBs as sources of UHECRs. As stated earlier, spectral ev- 
idence can distinguish between the two models, and, with 
the arrival of the next generation of detectors, the test 
is feasible (Sigl et al., 1995). Another significant observa- 
tional difference between the two models occurs when we 
consider the angular distribution of super-GZK events. In 
the coronal model, they must be isotropic because the GRB 
model (PRR) is designed to fit the observations (Meegan 
et al., 1992), except for the very small dipole anisotropies 
predicted, « a few times lO"'^. In the extragalctic model, 
super-GZK CRs must occur within a GZK sphere, with 
Rgzk ^ 30 Mpc. Within this distance the peculiar ve- 
locity of the Galaxy is formed (Scaramella, Vettolani and 
Zamorani 1994), so that we expect larger anisotropies. Here, 
I cannot help but notice the irony that it is the detection of 
an isotropy that would favor the local model, the only such 
case known to me in astronomy. 

For sufficiently large detectors, the total number of in- 
dependent directions of arrival in the coronal model must be 
at least as large as that of known GRBs, a few thousands. 
This differs sharply from the number expected for the extra- 
galactic model. Prom within a sphere R ^ 30 Mpc (corre- 
sponding to the path length of a CR which started out with 
E = 10^^ eV and is observed with E = 3x 10^" eV) we ex- 
pect in fact « 1 GRB every 300 years. The most reasonable 
time-delay over Rgzk is given by 

= 3 X lO^r f — 4— ^ ("t^V (4) 
" \^30 Mpc J \1 nGj ^ ' 

and this, assuming that the UHECRs coming from a GRB 

are spread over a timescale « M, means that ~ 100 inde- 
pendent sites from which super-GZK CRs are coming, are 
visible at all times. Thus about an order of magnitude sepa- 
rates the number of different acceleration sites visible in the 
two models, the coronal one being the more densely popu- 
lated. 

Assuming the validity of the extragalactic model, I 
would like to point out that the above-determined rate of 
GRBs inside the GZK-sphere, hazK = 1/10^ yr, allows a 
measurement of the time-delay between GRBs and UHE- 
CRs, in a different range than the one (~ lyr) discussed 
by Milgrom and Usov (1995a). In fact, suppose that experi- 
ments with perfect sky coverage were to reveal the existence 
of N small regions on the plane of the sky from which super- 
GZK events seem to arrive. Then, we can use the inverse of 
the rate above, l/hazK, as a universal clock to state that 
the time delay is « N/hazK- The areas of the small regions 
is limited from below by instrumental resolution, or by the 
deflections of particles along their flight-path to our detec- 



tors. Since, from conventional estimates of the magnetic field 
it is found that deflection angles of a few degrees are most 
likely for UHECRs (Sigl, Schramm and Battacharjee 1994), 
we expect to be able to resolve at most « 10^ such clusters. 
This means that the time-delay that can be measured is in 
the range l/hczK — '^0^ /ficzK, i.e., 10^ — 10® yr. A detailed 
study of this measurement will be presented elsewhere. 

In summary, the model of PRR, designed to fit the ob- 
servations of GRBs, also naturally accounts for the energy 
of the super-GZK Crs observed so far. The same model also 
avoids the efficiency problem of the extragalactic competi- 
tor, and is consistent with the statistical significance of the 
coincidences found by Milgrom and Usov (1995a). Angular 
distribution properties arc sufficiently distinct from those of 
the extragalctic model to make discrimination of the two 
models feasible with the next generation of detectors. 
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